Nonadiabatic transitions in electrostatically trapped ammonia molecules 
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Nonadiabatic transitions are known to be major loss channels for atoms in magnetic traps, but have thus 
far not been experimentally reported upon for trapped molecules. We have observed and quantified losses 
due to nonadiabatic transitions for three isotopologues of ammonia in electrostatic traps, by comparing the 
trapping times in traps with a zero and a non-zero electric field at the center. Nonadiabatic transitions are seen to 
dominate the overall loss rate even for the present samples that are at relatively high temperatures of 30 mK. It is 
anticipated that losses due to nonadiabatic transitions in electric fields are omnipresent in ongoing experiments 
on cold molecules. 

PACS numbers: 37.10.Pq, 37.10.Mn, 31.50.Gh 



The recent development of a large variety of methods and 
devices for the manipulation and trapping of neutral polar 
molecules offers new opportunities for molecular physics ex- 
periments Decelerated beams and trapped samples of 
polar molecules can be used to study intrinsic molecular 
properties, such as energy level splittings fS, ^ and life- 
times of metastable states ^ |5|, with unprecedented preci- 
sion. Once the densities of the trapped molecules become 
high enough and their temperatures become low enough, the 
inter-molecular interactions are anticipated to enable interest- 
ing new studies and applications mill- For all these studies it 
is not only of importance to increase the phase-space density 
of the trapped molecules but also to increase the time during 
which the molecules can stay confined in the trap, i.e., to re- 
duce the trap loss processes. Neutral polar molecules in low- 
field seeking states are routinely trapped in magnetostatic or 
electrostatic traps by exploiting the Zeeman or Stark effect, 
respectively |8jl9J. These traps typically exhibit zero field at 
the trap center. Within a certain area around the trap center, 
the trapped molecules can undergo nonadiabatic transitions, 
widely also referred to as spin flip or Majorana transitions, 
from a trapped state into a state in which the molecules are no 
longer trapped. In atomic physics, nonadiabatic transitions in 
quadrupole magnetic traps seriously hindered the generation 
of the first Bose-Einstein condensates, as these spin flips made 
it impossible to reach the required ultracold regime l lOllTll . 
The trap losses associated with the presence of the zero field 
at the trap center were eliminated by implementing the TOP 
(time-averaged, orbiting potential) trap on the one hand |12| 
and by keeping the atoms away from the trap center with an 
optical plug on the other hand |[T3| . Generally, trap loss due to 
nonadiabatic transitions can be completely suppressed by cre- 
ating a non-zero field minimum in the trap center. Akeady in 
1962, loffe introduced a special variation of a magnetostatic 
trap with a field offset in the center for nuclear physics exper- 
iments 1J4 |. Piitchard suggested in 1983 to use such a trap for 
the confinement of neutral atoms ifTsl . This type of magne- 
tostatic trap is now widely known as the loffe-Pritchard (IP) 
trap. 

For trapped polar molecules, losses due to nonadiabatic 



transitions have not been experimentally reported upon yet. 
Nevertheless, several possible geometries for an electrostatic 
analogue of an IP-type trap have been suggested to prevent 
these possible losses, such as a chain-linked trap lfT6l or a 
six-wire trap |17|. More recently, a trap with an electric field 
offset in the center has been demonstrated for Rydberg atoms 
ifTSl . Theoretical studies on the loss of molecules from purely 
electrostatic traps |19| or, more generally, from magneto- 
electrostatic traps 1 20 1 due to nonadiabatic transitions have 
also recently appeared. In both studies it was concluded that 
- although this can be different for any particular molecule - 
the loss rate is negligible at the milli-Kelvin temperatures that 
are currently achieved in electrostatic traps. 

In this Letter, we experimentally study trap losses due to 
nonadiabatic transitions in electrostatically trapped ammonia 
molecules. We quantify this trap loss mechanism from mea- 
surements of the trapping times in an IP-type electrostatic trap 
with either a zero field or an offset field at the center, under 
otherwise identical conditions. Moreover, we measured the 
trapping times for three isotopologues of ammonia, namely 
I'^NHa, i^NDg and i^NDg. These molecules have subtle dif- 
ferences in their energy level structure, resulting in different 
probabilities for nonadiabatic transitions. Contrary to the the- 
oretical expectations, the results demonstrate the enormous 
importance of trap losses due to nonadiabatic transitions even 
at the present 30 mK temperatures. 

Figure [T| a) shows the energy of the | J, K) — |1, 1) level of 
NH3 and ND3 in electric fields up to 150 kV/cm. The upper 
(lower) inversion doublet components of this level are seen to 
split into a low-field (high-field) seeking set of hyperfine lev- 
els labeled MK=-l {MK=+l) and into a set of hyperfine lev- 
els labeled MK=0 that are basically not shifted in the fields. 
The zero-field inversion splitting of ^^NHa is with Wi„i,=23.7 
GHz considerably larger than the Wi„i,=1.59 GHz of ^''NDs. 
At low electric fields, the Stark shift of the MK=-1 set of lev- 
els is proportional to E'^P^mv, and is therefore considerably 
less for ^''NHs than for ^''NDa in a given electric field with 
magnitude E. In Figure [TJb) the Stark shift of the upper in- 
version doublet component of ^'^NHa is shown for low electric 
fields on an expanded scale, such that the individual hyperfine 
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levels can be recognized. The quadrapole interaction of the 
^''N nucleus causes the observed main splitting in three groups 
of levels, labeled with the quantum number Fi ETl . The nu- 
clear spins of the hydrogen atoms cause small additional split- 
tings. In the electric field, the overall splitting behavior into 
(almost) unshifted levels belonging to MK=Q and levels that 
shift quadratically with E, belonging to MK=-l, is recog- 
nized. The individual hyperfine levels have a definite angular 
momentum projection quantum number Mp, where F is the 
total angular momentum. As the electric field scales linearly 
with the distance from the center in a quadrupole electrostatic 
trap, a distance scale (x, in ^m) appropriate for the trap used in 
our experiment is indicated on top of the Figure as well. In the 
lower two panels of Figure [T]the behavior of the hyperfine lev- 
els of ^""NDa and ^^NDa in electric fields is shown ||22l. Note 
that the electric field scale, and thereby the distance scale, is 
now a factor four smaller than for ^^NHa. For ^^NDa and 
^^NDa, the magnitude of the inversion splitting is quite simi- 
lar (1.59 GHz and 1.43 GHz, respectively) but their hyperfine 
structure is rather different. As the deuterium atoms have a 
nuclear spin of one, the hyperfine structure in the deuterated 
isotopologues of ammonia is generally quite complicated. In 
^''NDs, the hyperfine levels again split in three groups of lev- 
els due to the quadrupole coupling of the nucleus ll22l . 
The ^^N nucleus has no quadrupole moment, resulting in less 
(11, instead of 16 for ^"^NDa) hyperfine levels, but these are 
closer spaced and without a clear sub-structure for ^^NDa ID. 

Polar molecules are oriented by the local electric fields in- 
side the trap, and they normally remain in the same quantum 
state and adiabatically follow the field while moving through 
the trap. The adiabatic eigenstates can be quantized with re- 
spect to the axis along the electric field vector E and can thus 
be assigned the quantum number Mp. Nonadiabatic transi- 
tions can occur when the molecules can not follow the rapid 
change of the direction of the electric field when they pass 
with velocity u at a close distance x to the region of zero elec- 
tric field at the trap center When this rate of change, given 
by v/x, is larger than the energy difference between levels 
in the molecule, nonadiabatic transitions can occur. As v is 
in general not oriented along E, nonadiabatic transitions with 
AAIp=0 as well as with AMp = ±1 will be possible. In the 
case of ammonia, nonadiabatic transitions from the MK=-1 
to the MK=Q set of hyperfine levels will cause trap loss. The 
frequency of these transitions at a distance x from the trap cen- 
ter is proportional to x^ jWinv From this it is seen that nona- 
diabatic transitions can occur when the molecules come closer 
to the trap center than a typical distance h cx {Winv'^Y^'^ ■ 
Figure [T|b)-(d) curves of v jx are indicated (in MHz) for two 
different velocities (u=l, 10 m/s), relative to the highest fre- 
quency MK=^ hyperfine level. The crossing point of these 
curves with each of the MK=- \ hyperfine levels gives the cor- 
responding distance h for nonadiabatic losses from that partic- 
ular hyperfine level. Obviously, when there is a large non-zero 
field at the center of the trap, nonadiabatic transitions from 
the MK=-\ to the Af A'=0 levels will no longer be possible. 
Nonadiabatic transitions within the MK=-\ set of hyperfine 




50 100 150 200 400 600 

E (kV/cm) E (V/cm) 




50 100 150 50 100 150 



E (V/cm) E (V/cm) 

FIG. 1: (Color online) Energy of the \ J,K) — |1,1) level of 
^"NHs (dashed line) and ^"NDs (straight line) in fields up to 150 
kV/cm (a). In panel (b), (c) and (d) the Stark shift of the upper in- 
version doublet hyperfine components in low electric fields is shown 
on expanded scales for ^''NHs, ^^NDs and ^^NDs, respectively. For 
the calculation of the electric field dependence of the hyperfine com- 
ponents, see Ref. |22|. The v/x curves indicate the rate of change 
of the direction of the electric field for a molecule with velocity ?; (1 
and 10 m/s) passing the trap center at a distance x (upper axis). 



levels can still occur to some extent, although also these will 
be largely suppressed as only the magnitude of the electric 
field, but no longer the direction, changes rapidly when the 
molecules pass through the center of the trap. 

In the experiment (vide infra), we detect molecules in all 
low-field seeking MK=-\ hyperfine levels simultanuously. It 
is evident from inspection of the Stark curves shown in Figure 
[T]that each of these levels will have a different nonadiabatic 
transition rate. Therefore, when the trap loss is dominated 
by nonadiabatic transitions, a multi-exponential decay is ex- 
pected to be observed. The rate for nonadiabatic transitions 
depends on how often the molecule passes through the vol- 
ume with cross-section 6^ near the center of the trap. This 
rate can be approximated by fosc{b/ba)'^, where fosc is the 
oscillation frequency of the molecules in the trap and bo is the 
mean value of the impact parameter of the trajectories of the 
molecules. With a typical value for fosc of 1-2 kHz and with 
feo on the order of 0.2 mm, a value of b around 5 /im will lead 
to a loss rate due to nonadiabatic transitions on the order of 
1 Hz. 

The experimental setup, together with an expanded view of 
the IP-type electrostatic trap, is schematically shown in Fig- 
ure |2|a). A detailed description of the molecular beam ma- 
chine, and in particular of the deceleration of a beam of am- 
monia molecules, is given elsewhere |23|. Decelerated pack- 



3 



(a) 



detection laser 

decelerator 



valve skimmer 

hexapole 




(b) 




=FV„ 



(c) 

g ' 






— ND, 
--- NHj 








S 0.5 








|o.25 







^ ^ 1 1 "f" i 1 ' 1 ' 




I 

^0.75 
SO 




fl 0.5 




a 0.25 













-2 -1 





1 2 




distance (mm) 



FIG. 2: Scheme of the experimental setup, with an expanded view of 
the IP-type electrostatic trap. The voltages applied to the six cylin- 
drical electrodes are shown (b), together with the resulting potential 
energies for NH3 (dashed line) and ND3 (straight line) along a line in 
the horizontal symmetry plane for a non-zero (top) and a zero (bot- 
tom) electric field at the trap center (c). 



FIG. 3: (Color online) (a) Lifetime measurements for ^^NDs in the 
two trap configurations, both on a linear and on a logarithmic scale 
(inset). In the inset, straight lines are included to guide the eye. (b) 
Lifetime measurements for ^^NDa with original (Vjj /q) and reduced 
(Vj5/q/4) trap voltages on a logarithmic scale. The measurements 
have been performed with intervals of 20 ms and with 16 averages 
per timepoint. The non-zero field in the trap center is 16 kV/cm for 
all measurements shown. 



ets of either NH3, ND3, or ^^NDs molecules leave the 
decelerator in the low-field seeking MK=-l hyperfine levels 
of the I J, A') = 1 1, 1) level with a mean velocity of around 20 
m/s, and with a full width at half maximum (FWHM) veloc- 
ity spread of about 10 m/s. The decelerated packet contains 
appproximately 10^-10^ molecules and has a spatial extent of 
about 2 mm along the molecular beam axis and is about 2x2 
mm^ in the transverse direction at the exit of the decelerator 
1231 . The molecules are decelerated to a mean velocity around 
zero m/s upon entering the trap, which is centered 8 mm be- 
hind the decelerator. After a certain trapping time, the trap 
is turned off and the molecules are state-selectively detected 
using a laser based ionization detection scheme. 

The trap consists of six cylindrical electrodes, shown more 
clearly in Figure |2|b). To the four outer electrodes (3 mm di- 
ameter) positive and negative voltages (±Vq) are applied such 
as to generate a quadrupole field. The two electrodes (2 mm 
diameter) centered on the symmetry axis of the trap enable 
to create an additional dipole field. When voltages of oppo- 
site polarity (iVo) are applied to these electrodes, a non-zero 
electric field is generated at the center of the trap; by applying 
the same polarity (+Fd) to the center electrodes, a standard 
quadrupole trap with a zero field at the center is obtained. In 
Figure|2]^c) the potential energy for the Af K=-l levels of NH3 



and ND3 in the trap is shown along a line in the horizontal 
symmetry plane, both for a zero (bottom) and a non-zero field 
(top) at the trap center. In the latter case, the trap depth is 
actually limited by the occurence of four saddle points of the 
electric fields in between the electrodes. For Vq=10 kV and 
V/)=3 kV, the offset electric field at the center is 16 kV/cm 
and the trap depth is about 290 mK for ND3 and 130 mK for 
NH3. With no electric field at the center, the trap is consider- 
ably deeper, about 1.2 K for ND3 and 750 mK for NH3. 

Figure [3ja) shows measurements of the density of 
^^ND3 molecules in the electrostatic trap as a function of the 
time during which the trap has been switched on, both for 
zero and non-zero field at the center. Although the number 
density in the zero field trap is originally higher due to a bet- 
ter trap loading, it is also seen to decay significantly faster. In 
both cases, the trap is considerably deeper than the tempera- 
ture of the trapped ensemble of molecules; the latter has been 
determined from expansion measurements to be around 30 
mK. The molecules can therefore only leave the trap via col- 
lisions, via optical pumping due to blackbody radiation or via 
nonadiabatic transitions. Losses due to collisions can either 
originate from (in)elastic collisions of the trapped molecules 
with background gas or from cold collisions amongst the 
trapped molecules. At the present densities in the trap (10^- 
10^ cm^-^), collisions of the trapped molecules with each 
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other can be neglected. We have experimentally verified this 
by changing the number density of trapped molecules over 
one order of magnitude, without observing any change in the 
trapping time. The ^''NDs molecules in the low-field seek- 
ing MK=-1 hyperfine levels can be optically pumped to ex- 
cited ro-vibrational levels due to blackbody radiation from the 
room-temperature environment, leading to a calculated finite 
lifetime of about 7 s |24J. The purely single exponential de- 
cay of the number of ^"^NDs molecules in the trap with the 
non-zero field has a 1/e time constant of 1.9 s. This decay 
time results from the combination of optical pumping due 
to blackbody radiation and (in)elastic collisions with back- 
ground gas (background pressure 10^* mbar). By simply 
changing the polarity on one of the center electrodes of the 
trap, under otherwise identical conditions, a much faster and 
multi-exponential decay is observed (Figure [3]^a)); the addi- 
tional trap loss observed in the trap with zero electric field at 
the center is solely due to nonadiabatic transitions. From these 
measurements, it appears that the nonadiabatic transitions for 
^''NDs can be approximated with only two time constants. 
Referring to Figure [T|c), the fast decay can be attributed to 
nonadiabatic transitions within the Fi=2 set of hyperfine lev- 
els, whereas the slow decay can be attributed to nonadiabatic 
transitions from the Fi=l to the Fi=2 set of hyperfine lev- 
els. When the nonadiabatic transitions are modelled in this 
way (with the pre-factors for the two exponentials determined 
by the number of Mp levels), time constants of about 0.25 s 
and 2.0 s, respectively, are extracted for these processes. The 
errors on these 1/e time constants are less than 20 ms. 

We have performed a similar series of experiments for 
^''NHs (data not shown); this is actually the first time that 
this ammonia isotopologue has been trapped at all. In the zero 
field trap, a bi-exponential decay is also observed for ^''NHs, 
which can be explained in the same way as for ^"'NDa. After 
correction for the overall decay as measured for ^^NHa in the 
trap with the offset field, the time constants for the fast and 
the slow process are determined as 0.13 s and 3.5 s. The faster 
initial decay can be explained by the larger value of b whereas 
the slower decay from the Fi=l to the Fi=2 set of hyper- 
fine levels than for ^"^NDs can be rationalized by the slightly 
larger energetic separation and the reduced number of levels 

of "NH3. 

In Figure |3jb), the density of ^^NDs molecules in the trap 
is shown as a function of time. In this case, the decay in the 
zero field trap is seen to be almost single exponential, and the 
1/e time constant due to the nonadiabatic transitions is 1.2 s. 
A more pure single exponential behaviour can indeed be ex- 
pected as all low-field seeking hyperfine levels behave very 
similar in the electric field (see Figure [TJd)). When all volt- 
ages on the trap (and thus also the electric field offset in the 
non-zero field trap) are reduced by a factor of four, the life- 
time in the trap with the offset electric field stays the same 
(data not shown). In the zero field trap, however, the decay 
is not only faster due to the larger value of b but it is also 
multi-exponential as the grouping of the hyperfine levels is 
now relatively more important. 



As evidenced by the experimental results presented in this 
work, nonadiabatic transitions can be a dominant loss chan- 
nel in electrostatically trapped molecules, even at tempera- 
tures of tens of milli-Kelvins. These trap losses can be ef- 
fectively prevented by using an electrostatic trap with a non- 
zero electric field at the center, or, for instance, by orbiting the 
molecules off-center in a storage ring |25 1. The importance of 
trap losses due to nonadiabatic transitions depends on the de- 
tailed energy level structure of the molecule under considera- 
tion. Contrary to the situation in magnetic fields, there is not 
necessarily a degeneracy of low-field and high-field seeking 
levels for molecules at zero electric field. As a consequence, 
molecules in levels that are exclusively low-field seeking, like 
^"^NHs and ^''NDa molecules in the Fi=l hyperfine manifold 
of the |J,/-i:) = |1,1) level (see Figure [T]:b,c)), can be rather 
immune to nonadiabatic transitions. For some applications it 
might not be desirable to introduce additional offset fields. In 
that case it might be possible to select another isotopologue 
with a slightly different hyperfine structure to effectively sup- 
press the losses due to nonadiabatic transitions. 
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